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Abstract 

We consider the one-loop quantum corrections to the Nielsen-Olesen flux 
tube of finite length L, by imposing periodic boundary conditions. The 
calculations are based on a recent evaluation of these quantum corrections to 
the string tension of an infinite vortex. The finite length corrections are finite 
from the outset. If the computation is restricted to the zero modes we obtain 
the standard Liischer term 7r/3L for a closed string. The inclusion of the other 
fluctuation modes of Higgs and gauge fields, using the numerically computed 
trace of the Euclidian Green's function, leads to corrections that decrease 
exponentially with L. We present numerical results for these corrections, 
discuss their possible relevance, and the limitations of the approach. 
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1 Introduction 



The vortex solution of the Abehan Higgs model in (3 + 1) dimensions, known 
from superconductivity [1], has been introduced in particle physics by Nielsen 
and Olesen [2] as a possible model for strings. Indeed the authors show that 
the vortex can be related to the the bosonic Nambu-Goto string (see, e.g. [3]). 
This connection was mainly discussed on the classical level, the corrections 
due to quantum fluctuations were, there and later on, mostly considered 
already within string theory. Within the underlying quantum field theory 
the one- loop quantum corrections have been computed only recently ^ . The 
fact that the collective oscillations of the string are related to zero modes 
of the quantum fluctuations was already used qualitatively in Ref. P]. A 
detailed formulation of this connection has recently been presented in Ref. 
[5]; a further aspect has been addressed there: in quantum field theory the 
renormalization of the energy of collective fluctuations becomes part of the 
ordinary renormalization programme, renormalization requires the inclusion 
of the zero modes into the computation of the one-loop corrections to the 
string tension, they are not regularized and renormalized separately. 

While for this reason the contribution of the zero modes to the string 
tension of an infinitely long string cannot be quantified separately in quantum 
field theory, this is no longer so for the finite corrections which appear if one 
considers a string of finite length. It is the aim of this work to elucidate 
this aspect and to determine these corrections numerically. A vortex of finite 
length can either be constructed as an open string or a closed string. For an 
open string we would need to provide end caps, e.g., in the form of magnetic 
charges. It is hard to imagine how one could possibly compute quantum 
corrections to such a configuration. The same would hold for a closed string 
in the form of a torus. The technique developed in Ref. [4] can be applied, 
however, to a vortex of finite length with periodic boundary conditions. This 
can be considered as an approximation to a realistic closed string if the length 
of the string is much larger than its transverse extension. Such a computation 
will be the main subject of this article. 

The text is organized as follows: In Sec. |2]we recall some basic formulae 
of the model and its quantum fluctuations, referring mainly to Ref. [1] for 
all details; in Sec. |3] we develop the formalism for computing the quantum 
corrections to the energy of a finite string, i.e., those one- loop corrections 
that are not already included in the one-loop corrections to the string tension; 
explicit calculations, analytical and numerical, are presented in Sec. HI the 
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results are discussed in Sec. El we conclude with a brief summary in Sec. El 



2 The model 

The Abelian Higgs model in (3+1) dimensions is defined by the Lagrange 
density 

Here is the complex scalar Higgs field and 
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The particle spectrum consists of Higgs bosons of mass = 2Aw^ and 
vector bosons of mass = g'^v^. We denote the ration of these masses as 
^ = uih/Mw = 2A/(7^. The vortex solution [Tl |2] is defined, in the singular 
gauge, by the cylindrically symmetric ansatz □ 
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where r = y x"^ + y"^ and ip is the polar angle. The equations of motion for 
/(r) and A{r) can be solved numerically, see, e.g., Ref. 0. In terms of these 
functions the classical string tension takes the form 




dA{r) 
dr 



df{r) 
dr 



+ [A{r) + 1]' 



(2.6) 



The parameter dependence can be written in the form (Jd = {ti / g'^)h{^) , 
where /i(^) varies from 0.75 at ^ = 0.5 to 1.34 for ^ = 2. Here and elsewhere 
we use units such that mw = 1- 

The fiuctuations around this classical solution consist of those of the real 
and imaginary part of the Higgs field and those of the transversal, longi- 
tudinal and timelike components of the gauge field; furthermore the gauge 



We use Euclidean notation for the transverse components, so = = —Ai etc. 
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fixing introduces the corresponding Faddeev-Popov fields. In the fluctuation 
energy the Faddeev-Popov contributions cancel those of longitudinal and 
timelike components of the gauge field, so that these are irrelevant here [3]. 
The remaining fluctuations form a 4 x 4 coupled system. As the classical 
solution is independent of time and z the Euclidian fluctuation operator has 
the general form 

M,i = -{dl^dl)h,^Muo- (2.7) 
The transversal fluctuation operator A^_Ljj is identical to the one of the in- 
stanton in the 2-dimensional model; is has been presented in detail in Refs. 
[6l H]. The Green' s function of A^jj has, in momentum space, the formal 
representation 

C?i,-(x^, x^, k,iy) = }^ j,2^k^ + X2 (2-8) 

where and '?/'°(x_|_) are the eigenvalues and eigenfunctions of re- 
specively. The trace of this Green' s function has been computed numeri- 
cally in Ref. not via (12. Sp . but using a Jost function formalism adapted 
to coupled systems, see Ref. [7] and the Appendix of Ref. [8]. 

3 The energy of the vortex of finite length: 
basic relations 

We now establish the formalism for computing the quantum corrections to 
the energy of a vortex of finite length. As announced in Sec. [H we do 
this in an approximate way by imposing periodic boundary conditions in the 
longitudinal coordinate z. 

Starting point is the formal definition 

oo ^ 
n=—oo a 

where Ea and are the energies of the eigenmodes of the fluctuation oper- 
ators around the vortex and in the vacuum, respectively. The variable kn is 
the momentum in the longitudinal direction of the vortex, which for periodic 
boundary takes the values fc„ = 2tt/L. The energies of the eigenmodes have 
the form 

Eaikn) = VW+K (3.2) 
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and analogously for EQ^ikn). Here are the eigenvalues of the fluctuation 
operator A^^ in the transverse variables, as defined in Sec. HI 

We have introduced in the previous section the Green' s function of the 
fluctuation operator and its formal representation (12. 8p . We define a function 
F{kn, i^) as 

F{kn, ^) = j d^^^ Tr (^(x_L, x_L, fc„, v) - ^°(x_l, x_l, fcn, v)) , (3.3) 

Here Qq is the free Green' s function. We denote the eigenvalues of the free 
fluctuation operator by [Xa^'Y- We then obtain the relation 



(3.4) 



which is the basis of our numerical computation. In terms of F{kn, i^) we 
find for the energy of a string of finite length 

/oo J ^ 
F{27rn/L,u). (3.5) 

Actually Fikniv) only depends on fc^ + z/^. So with the definitions p = 
and F{p) = F{kn, v) we may write this as 



/OO J 



(3.6) 



with Pn = y ^"^ + (2nn/Ly. This can easily be converted into a weighted 
integral over F{p) via 

EniL) = - ri±w{p)Fip), (3.7) 

Jo 

with 

oo 

w{L,p)= J2 Vp' - {27in/LyQ{p^ - (2Tin/Lf) . (3.8) 

n=— oo 

The function F{p) has been computed numerically in Ref. [4j for various 
values of the parameter ^ = mH/mw- At small p it behaves as 2/p'^, due 
to the presence of the two translation zero modes. At large p it behaves as 
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a/p^, where the coefficient a is determined by the leading-order Feynman 
diagrams, see Eq. (7.1) of Ref. One easily convinces oneself that the 
weighted integral in Eq. (13. 7p is quadratically divergent. In the limit L — )■ cxo 
the sum over n can be replaced by an integral and one finds 

lim "^^^'^^ - ^ ^'^'-^ _ 1 



L—^OQ Zj 



1 1 TT 

— / dx^Jp^ -x^Q{p^ - x^) = — . (3.9) 
2tt J^^ 2n 2 



This limit yields that part of the string energy which is proportional to its 
length, E oc La and defines the string tension, or energy per length, of the 
vortex of infinite length ^ 

.. = -jr~£!^f(p). (3,10) 

Here we want to determine the corrections which arise for finite length. It 
is convenient, therefore, to subtract from Ef{{L) the term a^L . This can be 
done by redefining the weigth w{L,p) by subtracting the asymptotic weight: 

w,mp)=w{L,p)--p^= J2 Vp' - {27rn/Lye{p' - {27rn/L)^) - -p' . 

n=—oo 

(3.11) 

The finite length correction to the energy then becomes 

AE^{L) = E^{L) -a^L = - H l^w,{L,p)F{p) . (3.12) 

In this way we have gotten rid of the divergences. Indeed we have already 
used all counter terms of quantum field theory in order to obtain a finite string 
tension cTfi, as described in detail in Ref. [1]. One can verify explicitly, as we 
will see in the next section, that for a function F{p) which asymptotically 
behaves as c/p^ the weighted integral of Eq. f l3.12p is UV finite. The function 
Wsij)) does not tend to zero at large p, it oscillates with a period of Ap ~ 
27r/L, without being strictly periodic. The oscillations are due to the fact 
that more and more terms are included into the sum over n. Ws{p) is displayed 
in Fig. [H for L = 10. 



•^This expression for is of course divergent. F{p) has to be used in subtracted form, 
and the divergent parts have to be renormahzed, this is the subject of Ref. [4]. 
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Figure 1: The subtracted weight function Ws{p)- Sohd hne: Ws{p). The 
dashed hues indicate a square root behaviour of the maxima and minima. 



4 Numerical and analytical calculations 

As the function F{p) is know numerically from previous computations it 
seems a straightforward matter to evaluate the weighted integral of Eq. 
fl3.12p . However, this integral is subtle numerically, due to the oscillating 
and spiky weight function. It is useful, therefore, to begin with some related 
analytical and numerical calculations. 

The function F{p) is dominated, at small p by the two zero modes describ- 
ing collective oscillation, and it is instructive to compute their contribution 
to AEf[{L). For (Ta such a separate computation was not possible, as the 
contribution of the collective oscillations is infinite. When computing the 
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correction AEf[ this problem does not arise. We define the integral 



p dp 

Io{L,p) = I -^Ws{L,p) 



[Lp/27T] 



= 2 E 

+p- 
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vr L 
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^Jp^ — [2-nnl Ly — arccos 



Lp 



4 27r 

and its limit as p — t- oo 



P 



Io{L)= lim/o(L,p). 



In terms of Io{L) the contribution of the two zero modes is 

2 ^ 



A£;coii(^) = —h{L) . 

TT 



(4.1) 



(4.2) 



(4.3) 



At finite p the sum on the right hand side of Eq. (14. ip extends up to = 
[Lp/2T[\. It can be evaluated using the Euler-Maclaurin summation formula. 
One finds 



(4.4) 



We display the integral Io{L,p) in Fig. [2] for L = 10. One sees that the 
integral oscillates as a function of p, the width of these oscillations narrows 
only slowly, as l/^/p, and the band is slightly asymmetric. At p = 500 the 
width is of order 0.002 and the mean value can be read off, by taking the 
average of maxima and minima, with a precision of order 0.0001. For L = 10 
one finds Iq{L,p) ^ 0.1644; this is consistent to four digits with the analytic 
result 7r^/6L = 7r^/60 = 0.164493. The same value is found by evaluating 
the weighted integral via numerical integration. 
Using Eq. (14. 3 p we find 



AE,,n{L) = -71 /3L 



(4.5) 



This is the Liischer term [9l [TO] for a closed string (see, e.g., [H]). For a 
closed string the mode energies are 27m /L instead of vr/L for an open string, 
and the oscillating modes can propagate in two directions, up and down the 
z axis. This explains the factor 4 relative to the usually quoted value of 
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Figure 2: The integral Io{L,p) of Eq. (14. ip . Solid line: the function Iq{L,p); 
dashed horizontal line: the asymptotic limit 0.1644. 

7r/12L. We note that our result is obtained by evaluating an expression that 
is finite from the outset. This is analogous to finite temperature corrections 
which likewise do not need a regularization and renormalization. 

An important property of our weighted integral I{L,p) is apparent in Fig. 
121 the mean value between maxima and minima of the oscillations attains 
the final value of the integral already at very low values of p , i.e., after a few 
oscillations with period 2tt/L. This means in general, that in a region where 
F{j)) oc l/p^ the integral of Eq. f l3.12p will just oscillate around an almost 
constant average value. As F{p) ~ c/p^ for large p this means that in the 
asymptotic regime the average value of the integral reaches its asymptotic 
limit quickly. 

In order to perform the weighted integral for the realistic case of Eq. 
f l3.12p we need F{p) for a very narrow grid of values p, Ap -C '2,tt/L, as the 
weight varies strongly and the integration implies subtle cancellations. The 
numerical computation or Ref. |1] has provided values only on a relatively 
coarse grid of values p. As a computation of F{p) requires substantial CPU 
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time, a few minutes for each value of p, while the function itself is varying 
smoothly, it is convenient to use fits through the existing data points. One 
may use spline fits, but it turns out that a simple parameterization 



gives a surprisingly good global fit to the data, for all values of ^ = niH/mw- 
The ansatz can be understood as a fit where the effect of all higher modes, 
possible bound states and continuum, is simulated by just one pole with an 
effective degeneracy d at an energy A on the imaginary p axis, which contains 
the physical cut. The symbol A refers to Eq. fl2.8p . A^ being an eigenvalue of 
the transversal fluctuation operator. The parameters d and A^ 

for various values of ^, determined by a fit-by-eye, are given in Table [TJ 
The fits and the numerical data are displayed in Fig. [3] for three values of 
^ = rriH/mw. If necessary the approach could be improved systematically 
by including more poles into the ansatz fl4.6l) . This type of interpolation can 
be considered as a Fade approximation (see, e.g., Ref. [12j, Sec. 5.12) and is 
well adapted for functions with a cut in the complex plane. The parameters 
could be determined, e.g., by least squares methods. 

A considerable advantage of our fit is the fact that we can do the integral 
over p for the second pole term analytically, as well. We thus can avoid a 
very subtle numerical integration. One finds 



The delicate cancellations between the sum and the other terms on the right 
hand side can be avoided by subtracting Iq{L,p) and adding its asymptotic 




(4.6) 




n=l 




(4.7) 
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Figure 3: The fits to the function F{p) of Eq. (14. 6p . We display the data 
and the fits for p'^F{p) for ^ = 0.6 ,^ = 1.0 and ^ = 2.0. Symbols (+, x, *): 
data; solid lines: fits according to Eq. ( 14. 6p 



limit 7r^/6L of Eq. (14.41) . Of course the resulting expression 
h{L,p) = h{L,p)-Io{L,p) + 



7r2 



6L 

[Lp/2tt' 



2 E 



n=l 



2vrn 27m /- — — y(2Wl)M^ 

arccos — y {2Tin/Ly + A"' arccos 



L Lp p 



-Aarccos- + -A^ln^--— + — 4.8 
p y bL 

is not identical to Ix{L,p), but it has the same limit as p ^ oo. The oscilla- 
tions present in Ix{p, L) are considerably suppressed in I\{p, L) and the limit 
p — )■ oo, which we denote by I\{L), can be evaluated numerically without 
problems. To obtain the corrections to the Liischer term in Ai?fl(L), we have 
to multiply the result by —d/ir, so that, with our fit to F{p) we get 

^E,{L) = -^- -h{L) . (4.9) 
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In Fig. Hlwe plot Ix{L) as a function of L for = 0.7, 1.0 and 2.0, values 
which are in the range of the fit parameters given in Table [TJ The results 
obtained display a roughly exponential decrease with L. This is not unex- 
pected; had we imposed periodicity in time instead of periodicity in z we 
would have expected the thermal corrections to display such an exponential 
behaviour. We have to take into account that we do not consider the contri- 
bution of a single energy level, but the contribution of a cut generated by the 
pole in the transversal Green' s function. So the correction is not expected 
to be described by a simple exponential function aexp(— aL). However, the 
"effective" logarithmic slope a is in the expected range of values. For the 
range 2 < L < 6 we have a ~ 0.95 for = 0.7, a ~ 1.15 for A^ = 1 
and a ~ 1.55 for A^ = 2. The straight dashed lines in Fig. H] indicate this 
behaviour. We would naively have expected a = A = 0.837, 1.0 and 1.414, 
respectively. Indeed for L > 6 the slopes a decrease and may attain these 
values as L — )■ cxD. For L < 2 on the other hand the effective slopes increase. 




Figure 4: The integral I\{L). diamonds: A^ = 0.7; triangles: A^ = 1.0; 
circles: A^ = 2.0; the dashed lines indicate a simple exponential behaviour, 
as specified in the text. 
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5 Discussion of the results 



Using the results of the previous section the final results for AEa^L) can be 
easily obtained, using the parameters d and from Table [H We present, 
in Fig. [5] the function c(L) = LAEa{L) for ^ = 0.5 and ^ = 1.2. The 
function c(L) can be considered as the coefficient of an "effective" Liischer 
term c{L)/L. With this formulation we follow the presentation of lattice 
measurements of this term in QCD in Ref. |13]. As one sees from tabled] the 
values for d and A are 8.5 and 0.7, respectively, for ^ = 0.5, for = 1.2 they 
are 2.4 and 1.0, respectively. For the first parameter set the corrections to a 
pure Liischer term are sizeable even at L = 4 and have decreased to the 10% 
level at L = 6, while for ^ = 1.2 they have decreased to this level already at 
L = 3. For higher values of ^ the effective degeneracies become smaller and 
so does the deviation from c = —n/3. 
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Figure 5: The coefficient c(L). circles: ^ = 0.5; diamonds: C, = 1-2; triangles: 
A^ = 1.0; circles: A^ = 2.0; the dashed lines indicate a simple exponential 
behaviour, see text. 

There are various aspects under which we can consider the results of these 
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computations. It is satisfactory, at first, that the Liischer term appears in a 
straightforward way. As this term has been measured for QCD strings [13] 
it is certainly a part of estabhshed physics. The fact that the corrections 
due to higher fiuctuations appear here within the same formahsm and in an 
analogous way gives us confidence that these terms as well are not artefacts 
of the approximation but terms that would appear as a result of a suitable 
measurement. 

Of course our analysis is an approximation; we see two essential limi- 
tations: (i) we have to require L to be "much" larger than the transversal 
extension of the vortex and (ii) the corrections have to be small enough for 
the semiclassical approximation to be reliable. 

The transversal extension of the vortex is of the order of max( 1 / itih , 1 / my/). 
In our computation the mass my/ is set to unity, so the transversal extension 
is i? ~ max(l, 1/^); for the range of ^ considered here this is between 0.5 and 
2. An optimistic guess of "much larger" could he L > R and a pessimistic 
one L > lOR. In the second case our correction will be be negligible for all 
values of ^. The range of validity also depends on what is really measured. 
If, e.g., the corrections were measured on a lattice with the same periodic 
boundary conditions, then our analysis would be valid in the whole range of 
parameters considered here. The limitation of the approximation depends, 
unfortunately, mainly on the effects which we have neglected: the influence 
of end caps for an open string, and the effects of curvature for a closed vortex. 
Their magnitude will we difficult to estimate. 

The second limitation is the validity of the semiclassical approximation. 
The classical string tension is given by Tc/g"^ times a number close to unity. 
The one-loop corrections to the string tension were found, in Ref. [1] to be 
smaller than 0.5, so this correction is small even for g as large as unity. The 
corrections for a finite string, as found here, are much larger. Of course the 
ratio depends on L: ad is multiplied by L while c(L) is divided by L. But 
for g — 1 and L ~ 2 the ratio of these contributions to the string energy is 
not small. In any case, for L small enough the Liischer term and the even 
larger corrections will become comparable or exceed the classical energy and 
the semiclassical approximation will break down. For the lattice results it is 
found that the absolut e value of c{L) becomes sma//er than 7r/12 at small 
L, and possibly tends to zero. In Ref. [13] this behaviour is described by 
a relation derived using the QCD renormalization group and therefore relies 
on asymptotic freedom. Such an analysis does not apply here. The absolute 
value of our coefficient c{L) increases at small L and it is hard to see how 
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this could be different as the parameters d are positive throughout. A higher 
order resummation seems out of scope. 

6 Summary 

We have computed the one-loop corrections to the energy of a Nielsen- Olesen 
vortex of finite length. More precisely: we have computed the corrections 
that are not already included in the string tension. The latter were the 
subject of Ref. [1]. The corrections computed here are finite from the outset 
and, therefore, do not depend on renormalization conditions. The leading 
order correction at large L is the Liischer term which here takes the form 
7r/3L as appropriate for a closed string. The further corrections decrease 
exponentially at large L but can be, depending on the parameters, relevant 
for small L and intermediate L. Within our computational framework they 
appear on the same footing as the Liischer term and are related to the internal 
structure of the vortex. 

We have discussed briefly the limitations of the approach and conclude 
that, depending on the parameters of the model, there is a window in L where 
the corrections to the Liischer term are relevant and where their computation 
is rehable. 

We would finally like to point out that the method used here can be 
applied in a similar way to other vortex configurations, like cosmic strings 
(see [H] for a recent review). In fact the only information specific to the 
Abelian Higgs model was contained in F{p), the trace of the Green' s function 
of the fluctuation operator. Furthermore, the behaviour of F{p) at small 
p is determined by the zero modes and its asymptotic behaviour can be 
obtained from leading order Feynman graphs. So semi-quantitative estimates 
are easily accessible. 
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d 


A2 


0.5 


8.5 


0.7 


0.6 


6.4 


0.75 


0.7 


5.0 


0.75 


0.8 


4.1 


0.8 


0.9 


3.48 


0.8 


1.0 


3.0 


0.8 


1.1 


2.65 


0.9 


1.2 


2.4 


1.0 


1.3 


2.22 


1.1 


1.4 


2.06 


1.2 


1.5 


1.94 


1.3 


1.6 


1.87 


1.6 


1.7 


1.8 


1.85 


1.8 


1.75 


2.1 


1.9 


1.71 


2.35 


2.0 


1.68 


2.7 



Table 1: The parameters of a pole fit to F{p), Eq. f l4.6p 
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